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Abstract 

Our exact theory for continuous harmonic perturbation of a one dimen¬ 
sional model atom by parametric variations of its potential is generalized 
for the cases when a) the atom is exposed to short pulses of an external 
harmonic electric field and b) the forcing is represented by short bursts of 
different shape changing the strength of the binding potential. This work 
is motivated not only by the wide use of laser pulses for atomic ioniza¬ 
tion, but also by our earlier study of the same model which successfully 
described the ionization dynamics in all orders, i.e. the multi-photon pro¬ 
cesses, though being treated by the non-relativistic Schrodinger equation. 

In particular, it was shown that the bound atom cannot survive the exci¬ 
tation of its potential caused by any non-zero frequency and amplitude of 
the continuous harmonic forcing. Our present analysis found important 
laws of the atomic ionization by short pulses, in particular the efficiency 
of ionizing this model system and presumably real ones as well. 

PACS: 32.80.Fb, 03.65.Ge, 32.80.Rm, 02.30.-f 

We study a simple one-dimensional quantum system with the attractive 
potential, modeled by the (5-function. This system is assumed to be in the 
bound state until at some initial time t = 0 a) it becomes exposed to an external 
harmonic electric field or b) the strength of the system binding potential gets 
time dependent. The perturbation after a short interval T is turned off and 
our objective is to study the time evolution of the bound state on the interval 
0 < t < T. The pulses of external electric field are modeling the application of 
laser beams for atomic ionization. 

Excitation of the 5-function atom was studied in [1-6] for a simpler case of 
harmonic parametric perturbation of its potential when T = oo, i.e. of infinite 
duration, and the main conclusion was the complete ionization for arbitrary 
frequency and amplitude of perturbation. Other results have shown a surprising 
similarity of main features of the process with observed experimentally and 
numerically in spite of simplicity of this model. 

1. PROBLEM SET UP 

As in [1] we start by considering the one-dimensional stationary system lo¬ 
cated at x = 0 with an unperturbed Hamiltonian 


H 0 = - 


d*_ 
2m dx 2 


gS(x ), g > 0, —oo < x < oo, 


( 1 ) 


1 


which has a single bound state 


u b (p, s) = Vpe" p|11 , P=jp9- (2) 

In the continuous spectrum the eigenfunctions are 

u(p, k,x) = -±= (e ik * - , -oo < fc < oo, (3) 

V27T V P + #l ) 

with energies h 2 k 2 /2m while the bound state energy is W = —Hcoo = — h 2 p 2 /2m. 
Functions u(p , k, x), Ub(p, x) are normalized to 6(k—k') and to unity respectfully. 
Parameter m represents the mass of the bound charged particle. 

On the interval 0 < t < T acts a perturbing potential which is described by 
adding in Eq.(l) a time dependent term V(x,t) 

V(x,t) = eExrj(t), or V(x,t) = R5(x)i)(t), te[0,T], (4) 

where ri max = 1 and parameters E , R are responsible for the amplitude of 
perturbation. Thus we have to solve the time-dependent Schroedinger equation 

^ — Hoiptx, t) + V (a;, t)ip(x , t), t > 0. (5) 

at 

After expanding ip(x,t) in the complete set of functions u : 

/ OO 

0(fc, t)u(p, k, x)e~^ hk2/2m)t dk, t > 0, 

-OO 

( 6 ) 

the survival of the bound state at time t <T can be evaluated by |0(f)| 2 , if we 
assume the system to be initially in its bound state 

6>(0) = 1, 0(fc, 0) = 0. (7) 

It is more convenient [1] to proceed in dimensionless units (h = 2m = g/2 = 
1) and rewrite 


„«,(*) = m = - 1, »(*.*) = W 

where IT), is the rescaled energy of the bound state. The energies of states 
u(k, x) are W(k) = k 2 with multiplicity two for k ^ 0. These functions are 
normalized to 6(k — k'), while the bound state Ub(x) - to 1. 

2. DIPOLE FIELD PULSE PERTURBATION 

Beginning at t = 0 a perturbing potential Exrj(t) is applied to the atom 
and it stops at t = T. Here parameter E represents the electric field of the 
perturbation whose frequency is to. 

For solving Eq.(5) we use Eq.(6) and expand ip(x, t) on the interval (0, T) in 
the complete set (8) of functions u: 

/ OO 

Q{k,t)u(k,x)e~ lk t dk. (9) 

-OO 
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Then using their orthonormality and assuming cutoff of the perturbation po¬ 
tential for large |x| reduce dimensionless form of Eq.(5) to the following set 


m 


4Er](t) 



Q(k, t) 


e ~i(k 2 + l)t 

(fc 2 + l ) 2 


kdk, 


(10 a) 


G(k,t) = 


AEe(t)e i ^ 2+1 '> t sinojt k 
\[2 tx (k 2 + l) 2 


( 106 ) 


The evolution of Q(k,t ) is determined in Eq.(lOb) by 9(t) only. Using Eq.(7), 
then integrating Eq.(lOb) in time, and substituting the result into Eq.(lOa) we 
obtain a single equation which describes the evolution of 0(f) 


m = 


(4 Ef 

27r 


V(t) / 9(t')r)(t')dt' 


f-oo + 


(fc 2 +1 y 


-k 2 dk, 0 < t < T. 


( 11 ) 


The ionization probability of the bound state is 

P{t) = 1 - |0(f)| 2 , for t < T, 


( 12 ) 


and it becomes constant P(T) for all later times t > T. The internal integral 
over k in Eq.(ll) can be expressed in terms of Fresnel’s integrals 


S(t - t') 
where 



e ~i(k 2 + l)u 

(k 2 + l) 4 


k 2 dk 


3-4 u 2 
48 


A + iu 


3 + 


4 u 2 — 4 iu 


24 


B, 


A = 7r[l — 4>('\/7«)], B = e lu \/n/iu — A, u = t — i! > 0. 


(13) 


Numeric realization of ionizing by dipole pulses 

We solve Eq.(ll) for 9(t) by using the following technique, which is implied 
by a known method for the Volterra equations. The integral in time is approxi¬ 
mated by the summation over discrete equal subintervals defined by equidistant 
points on the interval [0,t]: 0,5,25, ...,N5 = t. As it is clear from Eqs.(7) and 
(11) 0(0) = 1, 0(0) = 0, then we have 9(5) « 0(0) + 50(0) = 1. In this way 
9(25) = 9(5) +59(5) where 9(5) is evaluated by the integral (11) without involv¬ 
ing the end point at t = 25. And so on by consequent computations of 9(j5) via 
the terms with numbers 0,1,2,..., j — 1. This technique realizes an approximate 
quadrature of Eq.(ll) and its precision becomes better when subintervals 5 get 
shorter. This approximation replaces Eq.(ll) with the sum 

0( n 5) = - { -^5 V (n6) 0(j5)r,(j5)S((n - j)5), (14) 

7=0 

which allows to find approximately 9((n + 1)5) = 9(n5) + 59(n5). We study in 
this work sometimes quite long pulses, they require for an acceptable precision 
long computing time as 5 should be very small. This can be helped by improving 
the transition from 9(n5) to 9((n + 1)5) by involving two derivatives of 0: 

5 2 

9((n + 1)5) = 9(n5) + 59(n5) + —9(n5). (15) 
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Differentiating Eq.(ll) we have 


m = 


m 

27r 


2 r 


r/(t) / 9{t')r)(t')dt' 


e -i(fc 2 +l)(t-t') 

{k 2 + l ) 4 


k 2 dk 


57T 

128 


d(t)r] 2 [t) — 777 (f) / 0 (f') sinwf'df' 


r°° i(fc 2 +i)(t—«') 


( fc 2 + 1) 3 


-k 2 dk 


By expressing the last integral over k in terms of Fresnel functions as 

Tr , , f°° e -*“( fe2 + 1 ) 9 1 - 4m + 4it 2 a 1 + 2m 

V(u)=j (fc 2 + i )3 = 7 r[l-$(Vm)]---+ e V™ ^—, 


Eq.(15) can be rewritten in the following form 


(16) 


o x2 Z7>2 • x n 

0((n+ 1)<5) = d(nd) H--— { — v{nS) ^ 9{jS)ij(j6)V((n - j)<5) 

^ i=i 


(17) 


57T 

256 


9(n6)if (nS) 


V{nd) + |^(n<5)] j)S)}, 

i=i 


convenient for numeric computation of the ionization probability. 
A. Sin-wave pulses sin(wf) 


The ionization of our atom by the sin-wave electric field excitation, which 
models the laser pulses, means that 77 (f) = sin(u;f). For illustration everywhere 
in this work the ’short’ harmonic pulses will have only five cycles, N = 5, this 
can be easily realized now by experimental techniques [7-10]. For u> < 1 even 5 
cycles of oscillations last a relatively long time, sometimes T > 100 and the use 
of 9(t) is necessary. The results are exhibited in Fig.l, where time f is measured 
in numbers of cycles of perturbation. 



FIG.l. Ionization probability P(f) caused by harmonic pulses 
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The plot structure shows that external frequency w is doubled in the process, 
which is suggested in some measure by Eq.(ll). When the electric field E = 3 
the complete ionization occurs practically after the third cycle of the pulse. 

All computations are done using Maple on the intervals T = 5 — with various 
ui and E , therefore T = 78.54,104.7,157.1 for u> = 0.4,0.3, 0.2 respectively. An 
acceptable precision requires S < 0.04, i.e. the sums in Eq.(17) will have up to 
4000 terms and the computations is time consuming as the sum is of recursive 
nature. 



Fig.2. Ionization probability by five cycled harmonic pulses with ca = 0.2 and 0.3 

Fig.2 confirms our observation that when oj is far from the resonance the total 
duration of the perturbation pulse can be more important than its frequency, 
experiments [7-9] confirm this effect in real systems. The same behavior is more 
visible in Fig.2b, where the electric field is smaller E = 0.5. 

For a very qualitative comparison with experimental results we note that in 
the case of Cesium atoms, whose external orbitals have radii about 2 Angstroms, 
where the electric field is about 30 GV/m. When we rescale \Wcs\ = 3.89 eV 
to \Wb\ = 1 this would correspond in Eq.(14) to E ~ 8 in the dimensionless 
units, therefore our E = 0.5 means roughly about 2 GV/m , which is a very high 
electric field in laser pulses, but reachable for present techniques. The case with 
the field strength E = 0.2 (in physical units this corresponds E ~ 0.8 GV/m 
for Cs or ~ 1.6 GV/m for Tungsten with |W| = 7.86 eV ) was computed too: 
P(T) ~ 0.053 for to = 0.2 and 0.027 for u = 0.3 while the shape of curves P(t) 
is similar to plots in Fig.2. These results are only for orientation and clearly 
should be considered as qualitative due to the limitations of our model. It looks 
that for smaller E the value of P(T) becomes proportional to ~ E 2 in agreement 
with Eq.(17). 

B. Pulsed dipole forcing 

Here we consider electric bursts acting on the model atom: rectangular 
r](t) = 1 and bell-shaped r](t) = 4 (t/T — t 2 /T 2 ) on the interval 0 < t < T. 
Though the bell-shaped pulse has the same amplitude as the rectangular one, 
but its ionization efficiency is much lower in Fig.3 because its total energy is 
smaller and more importantly it does not have high frequency harmonics. The 
ionization probability P(t) is an oscillating function but there is an important 
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difference between plots in Figs.3a and 3b. While the ionization probability 
by a rectangular pulse of duration t\ is given by P{t\) in Fig3a, the ionization 
evolution for bell-shaped pulses is presented by Fig.3b, but only for the pulse 
of length 10, i.e. if T = t\ function P(t\) in Fig.3b does not represent the final 
probability and the corresponding computation should be performed namely for 
T = t\ because <5 = S(T). Our calculation with E = 0.2 and the bell-slraped 
pulse is not plotted in Fig.3b because the curve P(t) runs very low: its maximum 
P(4) ~ 0.013 and P(10) is less than 0.002. 



FIG.3a. Rectangular pulse FIG.3b. Bell-shaped pulse 

Fig.3. Ionization probability caused by rectangular and bell-shaped electric bursts 

Our approach is modified below for studying the ionization caused by the 
parametric modulations of the binding potential. Though this is hardly achiev¬ 
able in practice but, as we already mentioned, it exhibits some illuminating 
features of the ionization by short pulses which are quite universal. 

3. PULSED MODULATION OF BINDING POTENTIAL. 

Here we consider short pulse bursts of the potential strength which has to 
be studied by a different computations technique. Eq.(5) now has the following 
form 

= Hoil)( x ,t)+RS(x)ri(t)ip(x,t), r](t) = 0 when t 0 [0,T], t> 0, (18) 

with the initial conditions given by Eq.(7). Using the expansion (6) in terms of 
the stationary eigen-functions and the methods developed in [1] yields a simple 
equation for the function 

6(t) = 1 + 2 i [ Y(t')dt ', (19) 

Jo 

where Y(t) is defined by the following integral equation 

Y(t) = Rrj{t) jl + [2i +M(t-t')]Y(t')dt'^ . (20) 

The function M(s) in (20), see [1], is 

, . 2 i f°° u 2 e~ ls ( 1+u ) . / i _ is . p-. 

M(s) = — / -r- du = —i+\ —e +?.$(vfs). (21) 

7T Jq 1 + tr V 7TS 
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It behaves as \Ji/ns — i when s —> 0 and is proportional to s 3 / 2 e ls when 
s —> oo. 

A. Case of rectangular pulse 77 (f) = 1. 

For approximate evaluation of the function 9(t) in this set up we have to 
solve numerically the integral equation (20) that can be done if there is an 
effective way of computing M (s) for not very large values of s as our pulses are 
not long. Using [11] equation for the integral in Eq.(21) can be written in terms 
of gamma functions 



sfxe~ isx 
1 + x 


dx = e is T 



r 



and then one can apply the power series expansion for the confluent hyper¬ 
geometric function and present M (s) in a rapidly convergent form 


M(s) 


—i - 


l~i~ (— is) n 

V S7T 2 —/ (2 n — 1 ) n\ 

n =0 v ' 


( 22 ) 


For a given pulse duration T the upper limit in the sum (22) can be only slightly 
larger than eT to give a good precision for M(s). 

Our next step is express this function by a power series using the behavior 
of Y(t) near zero, 

Y(t) = ]T Cmt m/2 , 0 <t<T, (23) 

m =0 

and substitute Eqs.(22) and (23) into the integral equation (20). The result 
reads 


i?_1 Cmtm/2 = 1 +i J 2 ( 

m —0 m=0 


^l+m/2 

1 + m/2 


E B 

k,n—0 


k , 1 

2 + l’" + 2 


CkQ'nt 


n+- 


where B(x,y) = T(x)T(y)/T(x + y) is the /3-function [11]. Thus we obtain 
equations for finding coefficients c*: 


Cm 


c 0 = R, 



Raoco = 2 R 2 



2 iR 


L^J 

Cm -2 + R ^ B 
m 

71=0 


m 


+1 1 

2“- n ’ n+ 2 


^n^m- 2n— 1 1 


(24) 


m > 2. 


The symbol |_ccj as usual means the integer part of x, and coefficients c m in (24) 
can be calculated consequently as they are expressed in terms of c* evaluated 
earlier. Clearly we will keep in (22) and (23) only finite numbers of terms. 

By combining Eqs.(19), (23), and (12) we find the ionization probability as 
a function of the pulse duration T 


P(T) = 1 - 1 + 4i V -£l!—T 1+m/2 
K ^m +2 


771=0 


(25) 


The results of computation by Eqs.(22-25), where we keep about 80 terms in 
sums, are presented in Figs.4 which show this probability for different values 
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of pulse amplitude R and T < 10. Left Fig.4 shows also that a very short 
attractive pulse with R = —0.5, T ss 0.9 is more effective for ionization than 
longer pulses of the same amplitude. 



FIG.4 Plots of ionization probability P(T) for rectangular pulses 

In agreement with the common sense when the external pulse decreases the 
binding energy, R > 0 in (4), the ionization rate is larger. There is always an 
intermediate pulse duration corresponding to maximum of ionization and for 
longer positive pulses the ionization rate approaches to a constant. Right Fig.4 
exhibits the ionization probability caused by the attractive positive pulse whose 
amplitude is one half of the binding energy and duration T < 4.5. One can see 
the resonances and quite a high probability 0.32 when T = 0.55 because the 
pulse fronts have infinite slopes and thus introduce very high harmonics though 
the total pulse energy is limited. Note also that R should be compared with 
g = 2 in our units, |i?| = 1 is the largest considered here. 


B. Bell-shaped pulse rj(t) = A(t/T — t 2 /T 2 ) 

On the interval 0 < t < T we consider a pulse, symmetric about t = T/2, 
whose maximum amplitude is R , see Eq.(4). In this case the same approach 
as before converges slower and one needs more terms in Eqs.(22) and (23), we 
used up to 400-500 of them. 

A straightforward analysis of Eq.(20) shows that the power series expansion 
for Y(t) has the following form 


Y(t) = Y c m t 1+m/2 , 0 < t < T. (26) 

m —0 

By substituting (26) into Eq.(20) and using Eq.(22) we solve Eq.(26) and find 
the coefficients which define Y ( t ) 


Tc 0 
4 R 


1, ci = 0, 


TC2 
4 R 


1 Tea 
T' 4 R 


o-o cq B 



TC4 

~4R 


i—+a 0 c 1 B 



^ = i^jr- + a 0 c 2 B Q, 3^j + a-xCoB Q, 2^ - ^ a 0 c 0 B Q, 2^j , 


( 27 ) 



_ 2ic m _4 2ic m _g / 1 777. + 1 \ 

~±R~ “ T(m - 2) a ° Cm ~ 3 y 2 ’ 2 J 


+ E 


Cm—2n—5 


n—0 


2n ±1 


-«n+l ~ ^7 







?n > 6. 


One can see again that as in Eqs.(24) each coefficient c* is defined via corre¬ 
sponding ones found earlier. 



Ionization probability by bell-shaped Time evolution of ionization 

pulses of different duration T when T = 2.5 and R = ±0.5 


Note that in the case of a rectangular pulse after solving Eqs.(24) the solution 
can by found by using Eqs.(19) and (20) for any pulse duration t = T or in 
Fig.5a. But for the present bell-shaped pulse the meaning of plots in Fig.5b is 
different, see Part 2B, and also the ionization probability is found using Eqs.(26), 
(27) and (25) in the slightly modified form 


P(T) 


1-\6(T)\ 2 = 1- l±4i^ 


m =0 


C m rp2+m/2 

m ± 4 


2 


(28) 


where all Ck , k = 0,1,... depend on T. 

P(T) in Fig.5a is the ionization probability at the end of the corresponding 
pulse with R = —0.5 and length T, but the plot for maximum Max(P) which 
occurs at an intermediate time 0 < t < T and for longer pulses it is much 
larger than P(T). Fig.5b shows P(t) behavior in the case T = 2.5 when the 
pulse amplitudes are R = ±0.5. Positive pulses are more effective as before: 
P(T) = 0.13 is close to Mcix(P) when T = 2.5 while P{T) and Max(P) are 
0.139 and 0.147 respectively for T = 3. This probably means that change of 
sign of the external force plays an important role and might suggest that more 
realistic perturbation (a harmonic one with some envelope) can be more efficient. 
One can see that the system can be ionized before the end of perturbation pulse, 
i.e. the value of Max(P) is interesting. 

C. Short sin-wave pulses ?y(t) = sin(wt) 
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Our pulse exists on the interval 0 < t < T and as in Part 2 its length 
has always an integer number N = 5 of cycles, T = 2nN/u>. We assume 
in computations that the frequency of principal harmonic of Y (t) is u> and 
choose an integer K of harmonics sufficient for modeling Y(t). Then using the 
Galerkin [12] method for solving Eq.(20) the function Y(t) is approximated by 
the following sum 


K 

Y(t)= ^2 a kfk(t ), where f k (t) = e lujkt . (29) 

k=—K 


The solution method requires the discrepancy of using the approximation (29) 
in Eq.(20) be orthogonal to all functions f) c ■ This procedure creates the linear 
algebraic system for coefficients a*. 


£>* / dtf m (t)lf k {t) - Rr](t) f [2i +M(t')]f k (t-t')dt'\ = 


k=—K 


(30) 


R / rj(t)fm(i)dt , —K < m < K. 


By substituting Eq.(29) into (30) this system can be rewritten in the stan¬ 
dard form 

K 

^ v — B m , (31) 

k=-K 

where all coefficients with k ^ in ± 1 defined by the following relations: 


r _ . R_ ( M-k — Mi_ m M_fc — M_i_ m ^ 

k,m k-m 2w\^ k — m + l k — m— 1 ) 


(32) 


For k = m — 1 and k = in + 1 we have respectively 


R 

Cm-i,m = — [(1 + 2ituT)Mi- rn — M-i-m — 

4 to 

R 

C m +i,m = — [(1 ~ 2 iu>T)M_i_ m — Mi_ m + 2iu)M_i_ m ]. 
4 uj 


(33) 


Other terms in Eqs.(31-33) are 
0, k ^ m, 


Ak—m 


T, k = 


R 

B m — ^rr[-^l—m 

2i 


(34) 


M n = [2 i + M(t)}e inut dt, M n = / [2 i + M(t)]te inut dt. 


Our computation by Eqs.(31-34) for the cases when T = 5, R = 0.2, 0.4, 
and the pulse consisting of 5 cycles (therefore w = 2nN/T ss 6.3 > wo) are 
shown in Fig.6: 
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FIG. 6 . Ionization probability P(T) caused by sin-wave pulse 

One can see that with pulse amplitude 0.4 the ionization is quite effective. We 
used only five harmonics, I\ = 5, in computations but checked the precision by 
taking K = 10 which produced the curves P(t) almost identical to ones in Fig.4. 
The roughly linear time dependence of the ionization probability is caused by 
the fact that ui is significantly larger than uiq, i.e. the energy of ionizing photons 
in Fig.4 exceeds the binding energy of the model atom and thus we observe the 
first order effect. In reality this process might be compared with the soft X-ray 
ionization. 

For solving Eq.(20) this method of approximation Y(t) by Eq.(29) and the 
routine defined by Eqs.(31-34) becomes inefficient when oj is much smaller be¬ 
cause the interval T is too long in this case and oscillations of M(s ) are not well 
modeled by harmonics of 77 (f). We use the special properties of the Volterra 
Eq.(20) to model function Y (t) by a set of its discrete points. When the ioniza¬ 
tion is provided by at least two photons and oj = 0.6, 5 cycles of sin wave make 
T > 60, and thus one needs about N = 10 3 points for a decent approximation 
of Y(t). Denoting temporarily F(s) =2 i + M(s) the integral equation (20) is 
replaced by the following one, where we keep the integral term in the interval 
of divergent behavior of M(s) 


R Y (t n ) = r)(t n ) + rj(t n ) 


n—l ~A 

A ^ F(t n - t m )Y(t m ) + / F(s)Y(t n — s)ds 

m=l Jo 


(35) 

Here t n = nA, n runs from 1 to N and A = T/N. The integral in Eq.(35) is 
approximated using Eq.(22) for M(s) when s is small and only the linear term 
of Y (t)-dependence on the interval (t n _i,f n ): 


F(s) 


s 2 «3 

1 + ZS + y-*30 


Y(t n - s) « Y{t n ) + 


Y(t n -i)-Y(t n ) 


(36) 

Near small s we keep more terms in singular F(s) than in Y(t n — s) because Y(t) 
has a regular behavior there. By evaluating the integral in Eq.(35) we come to 
the recurrent equation for computing the set Y(t„) consequently starting from 
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n = 2. 


Y{t n ) 


Rv(t n ) 

l-B-iA/2 



n —1 

1 + A y ' F(tn — tm)Y(tm) + (A + iA/2)Y(t n -i) , 

m=l 

(37) 

A 2 *A 3 \ [iA / 4 4iA 2A 2 2*A 3 \ 

2T'l35 7’ “ V V ^3 + "lfT + T05 ~~ 945 ) ' 


Here we will neglect terms of the order A 6 and higher. 

It is clear that Y(0) = 0 and for approximating Y(t\) by a polynomial with 
the same precision we substitute F(A — s) into Eq.(20) and using Eqs.(36,37) 
obtain the following relation 


10 / 

^2 c k A 1+k/2 = ElwA- 


w 5 A 5 \ 

J 



iA2+n/2 

2 n/2 



where H n ( A) = 


B 


4 + n 1 
2 ’ 2 




- n 3 
_ ’ 2 


A 2 „ 

+ 6 S V 2 


4 + n 5 
'2 


(38) 

, 7 -) 

30 V 2 2 ) 


These eciuations are sufficient to evaluate c k , k = 0,10 and find Y(t\) = 
Y (A). We applied Eqs.(35-38) to compute the ionization by the five cycle pulses 
of lower than in Fig.6 frequency ui = 0.6 and amplitudes R = 0.2, R = 0.4. 



FIG.7. Ionization initiated by sin-wave pulse with u> = 0.6 

The dynamics of the process is shown in Fig.7. One can see that the ionization 
probability is higher than in the case of u> ss 6 in Fig.6. This result is probably 
caused by a much longer (more than in order) action of the perturbation while 
its average force (determined by R) roughly the same. The sin-wave pulses 
appeared to be quite efficient for ionization and P(T) seems to be proportional to 
R 2 for smaller R like for the laser pulse perturbation. The time of perturbation 
in Fig.5 is measured in the number of harmonic cycles, the pulse ends when 
N = 5 and T « 52.4. 

In the dimensionless units when the atomic ionization energy is 1 our results 
can be tentatively mapped onto real system again, say Tungsten W and Cesium 


12 





Cs with binding energies 7.86 eV and 3.89 eV respectively. It is easy to see that 
lo = 0.6 would correspond the wavelength about 260 nm for W and 530 nm 
for Cs in our cases. In experiments often are used ~ 4 — 10 fsec laser pulses of 
A = 800 — 830 nm, [7-10], and having in mind a qualitative application of our 
theory we perform a somewhat less precise computation for smaller lo = 0.4 and 
0.2 respectively when the pulse time T is longer. 



FIG.8. Ionization by short pulse harmonic waves of lo = 0.2 

The results presented in Figs.8 and 9 confirm the importance of total pulse 
duration. As before our pulses have only 5 cycles, the time of their action is 
T =157 dimensionless units when lo = 0.2 in Fig.6 and twice shorter when 
u) = 0.4 in Fig.9. Fig.8 roughly imitates the Cesium ionization when R = 0.1 
and 0.2 while A ~ 800 — 830 nm. The amplitude R = 0.2 is sufficient for almost 
complete ionization even at the pulse beginning. 

Fig.9 gives a hint of the Tungsten ionization (keeping in mind the same 
A ~ 800 nm). 



FIG.9. Short pulse ionization when atomic binding energy is larger 
The ionization level is lower than in case of Cs though the wave length is twice 


13 



shorter. This clearly agrees with greater binding energy in W atoms. 

As the parametric perturbation acts directly on the binding energy it is 
more efficient in Figs.7-9 than the more realistic perturbation by the external 
harmonic electric field in Part 2. These results show that the ionization of our 
model atom in some measure describes qualitative behavior of real systems. 

4. SUMMARY 

The atomic ionization by short pulses of external forces is studied on a sim¬ 
ple one-dimensional model which allows to construct an exact theory of the 
process and realize its conclusions by several methods of numerical computa¬ 
tions. This creates a basis for comparison with approximate solutions of more 
realistic models, simulations, and experiments. Our main results include the 
observation that for external frequencies, much lower than the resonance ones, 
the total duration of the pulse is more important for effective ionization than its 
frequency. When the ionization level is substantially far from the complete one 
it is increasing approximately linearly in time and has resonances as a function 
of pulse duration. For ionization caused by the dipole electric field the frequency 
of these resonances is twice larger than the frequency of external forcing. 
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